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Is it possible to prove that the properties of Hamiltonians, such as the ground state energy, results 
of dynamical evolution, or thermal state expectation values, can be efficiently calculated when the 
Hamiltonians have physically motivated constraints such as translational or rotational invariance? 
We report that rotational invariance does not reduce the difficulty of finding the ground state energy 
of the system. Crucially, the construction it preserves the translational invariance of a Hamiltonian. 
The failure of the construction for the properties of thermal states at finite temperatures is discussed. 



Introduction: In order to study the quantum systems 
that appear in nature, it is desirable to simulate their 
properties on a computer. A number of techniques, such 
as Density Matrix Renormalization Group and Quantum 
Monte Carlo simulations appear to work in a wide range 
of interesting cases for calculating ground state energies, 
dynamical evolutions and expectation values of thermal 
states. However, none of these schemes are guaranteed 
to perform efficiently. Can such guarantees be made? 
Recent results in the field of quantum information have 
allowed a classification of the difficulty of some of these 
problems. The dynamics of a ID Hamiltonian can im- 
plement any quantum computation pQ, and thus there 
must be instances which are hard to simulate on a classi- 
cal computer (BQP^P). Finding the ground state energy 
of a one-dimensional Hamiltonian with nearest-neighbor 
interactions between 12-level systems has been shown to 
be QMA-complete [2J, which means that the result can 
be checked to be a solution only on a quantum com- 
puter, an even harder task (QMA^BQP). This result 
represents the culmination of a sequence of incremental 
steps including [31 IH [5] . Under the restriction that the 
energy gap between the ground and excited states is con- 
stant for increasing system size, expectation values can 
be efficiently evaluated classically [TJ, meaning that this 
restricted problem is within NP (where a solution can be 
verified on a classical computer). 

One might ask whether other restrictions, especially 
those that are physically motivated, make the ground 
state energy, or thermal state expectation values, easier 
to find. Those that are of particular interest, as applied 
to an iV-spin Hamiltonian H , are translational invariance 
(TI) and rotational invariance (RI), 

TH1* = H; U® N HU^ N = H 

where T is the permutation operator T \i\i2 ■ ■ ■ ijss) = 
1*2*3 ■ • ■ ijv*i)i each i runs over the full dimension of the 
local Hilbert space, and the RI condition is required to 
hold for all single spin unitaries U. Perhaps such restric- 
tions would allow us to guarantee that classical simula- 
tion is possible. For example, RI quantum states have 
two-body reduced density matrices which are Werner 
states. Given that these are a function of a single param- 
eter, this indicates a reduction in the parameter space, 



thereby motivating such hopes. Nevertheless, it has 
been shown that the dynamics of a TI and RI Hamil- 
tonian remain BQP-complete [HUH] and, similarly, QMA- 
completeness can be retained in a translationally invari- 
ant Hamiltonian if the interactions are allowed to extend 
over a logarithmic number of spins |10j or non-locally 
[8]. However, for estimating the ground state energy of 
a local TI Hamiltonian, complexity theoretic arguments 
are not applicable in the same way as the non-TI case. 
The discussion of the complexity of a function f(x) is 
based on how difficult it is to evaluate f(x) as the size 
of the input increases. In particular, if the time required 
is upper bounded by a polynomial of [log x] = N, then 
the problem is said to be in P i.e. it requires polyno- 
mial time. Roughly speaking, in the case of ground state 
problems, the nature of the Hamiltonian terms replaces 
the function input, and the size of that input is related 
to the number of qubits N (for example, for a nearest- 
neighbor Hamiltonian, one should specify a term for each 
neighboring pair of qubits). However, in the case of a TI 
system, there is no parameter that scales, except for N 
(which is of size log N, so solutions that are polynomial 
in N are in the class EXP). Instead, for TI problems, we 
are interested in how long the algorithm takes to run as 
a function of N. It has recently been proven that this 
running time can only be polynomial if two complexity 
classes, BQEXP and QMAexp are equal [11 J. This was 
achieved for a nearest-neighbor Hamiltonian with a ID 
geometry, although the local Hilbert space dimension of 
the spins is extremely large. 

In this paper, we explore the influence of RI on ground 
state and thermal state problems, and show that for 
ground state properties, the case of translational and ro- 
tational invariance (TRI) is of equal difficulty to that of 
TI, i.e. rotational invariance does not reduce the com- 
plexity of ground state calculations. 

Decoherence-Free Subsystems: Central to any discus- 
sion of rotational invariance is the subject of decoherence- 
free subspaces and subsystems [TH [T31 dH [IS], which 
allow one to specify how to encode information, cither 
quantum or classical, against globally applied noise of 
the form U® N , for some allowed set of unitaries U. An 
excellent review on the topic can be found in |16j . For 
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concreteness, consider encoding information on qubits, 
taking U to be any SU(2) rotation; similar construc- 
tions can be achieved for d-dimensional systems subject 
to any SU(n < d) rotational symmetry, or discrete subset 
thereof [T71 ES]. 

The 2 r -dimensional Hilbert space of r qubits can be 
decomposed into the form 



r/2 



where the sizes of the subsystems are \A4j \ — 2j + 1 and 
\J\fj\ = £ for even r, and \J\fA = c r r\ + c r ~\ for odd r 

(note that j is an integer for even r, and a half-integer 
for odd r). The function is defined by 



2j + l 



r/2-jJ r/2 + j + l' 



The subsystems Afj are referred to as 'decoherence-free' 
i.e. they are not affected by collective rotations. Classi- 
cal data can be safely encoded across all such subsystems. 
Since coherence is not preserved in superpositions across 
different subsystems, quantum information must be en- 
coded in a single subsystem. In the limit of large r, the 
largest of these is for j — y/r/2, with size 2 r ~ log 2 r . P r ^ 
will be used to denote a projector of r qubits onto the 
j\fj subsystem, reserving subscripts for denoting which 
qubits the projector is to be applied to. 

For any to > r/2+j qubits that are a subset of r qubits 



r,j nm,m/2 



pr,3 p 



= 0, 



(1) 



This comes from the addition of angular momenta - the 
block of to spins has total angular momentum to/2, and 
the remaining spins cannot have more than (r — to) /2 an- 
gular momentum. It is impossible to get a total angular 
momentum smaller than the difference, which is neces- 
sarily greater than j if to > r/2 + j, and hence the two 
projectors have overlap. By a similar argument, 



pm,m/2p2,Q q 



(2) 



The decompositions for the 3- and 4-qubit cases are given 
in |17j . and one can readily verify, for example, that 
p4,0p3,3/2 = o or P 4 ' 2 P 20 = 0. 

Hamiltonian Properties: Given the ability to encode 
information in decoherence-free subsystems, which have 
dimensions approaching those of entire quantum systems, 
it is relatively simple to convert Hamiltonians into RI 
Hamiltonians. Evidently, the key is to preserve the im- 
portant features of the initial Hamiltonian such as which 
state is the ground state, the energy gap, locality of the 
Hamiltonian, other symmetries etc. 

Lemma 1. A Hamiltonian H\ — hi of N spins of lo- 
cal Hilbert space dimension d interacts with k-body terms 
which are local on a ID lattice, where < 1, d and k 



are 0(1), and the gap between the ground state and the 
first excited state is A. There exists another Hamiltonian 
H 2 = J2i K °i 'O(N) qubits interacting with k 1 -body local 
interactions and gap A that has the same ground state 
energy and is RI, with k' ', \h[ \ ~ 0(1). 

This problem does not, yet, incorporate an assumption 
of TI. With no restrictions on the energy gap A, it is 
known that the ground state problem is QMA-complete 
for fc > 2 0, and this also holds when the lattice is 
restricted to a spatial geometry, such as ID (2j |9]. On 
the other hand, if A is a constant, then Hastings has 
proved that the problem is within NP [7J. The following 
proof is only restricted to the ID case for pedagogical 
reasons, but the lemma applies to arbitrary geometries. 

Proof. We transform any instance of a Hamiltonian Hi 
into a RI version H 2 with the same energy gap A by 
finding a decoherence-free subsystem of r qubits of size 
d r > d. We now assume that Hi was composed of 
spins of dimension 2 r rather than d. To ensure that 
the ground state and energy gap are unaffected by this 
mapping, we apply a local magnetic field J^^^ 1 N) 
J = 2/cmaxi \hi \ to every site to shift all other states to a 
higher energy above the ground state. Having done this, 
one can trivially map the levels to d r — 1 of the spin to 
the d r levels in the decoherence-free subsystem. Hence, 
Hi acts on rN qubits with k' = rfc-body terms which are 
local on the same lattice, and the energy gap is the same. 
Our reason for selecting J in such a way is that were we 
to take the original ground state and remove one of the 
spins, the total energy cannot be more than fcmaxj \hi\ 
below the ground state energy, so replacing it with a state 
from one of the other subspaces must increase the energy 
over the ground state by at least fcmax^ \hi\ > A. □ 

It is important to observe that by encoding d- 
dimensional spins in the decoherence-free subsystems of 
qubits, expanding the number of spins we interact with 
(k' > k), we lose the ability to reassert nearest-neighbour 
interactions by recoding into a larger local Hilbert space 
dimension - this would break the rotational invariance. 

For example, the difficulty of finding the ground state 
energy of a RI Hamiltonian with 16-body interactions is 
the same as that of a non-RI Hamiltonian with 12-level 
systems and nearest-neighbor interactions, i.e. QMA- 
complete [2]. Given that we are interested in classical 
simulation, we must ensure that the encoding process of 
Hi — * i?2 can be efficiently simulated. This is guaranteed 
because the encoding process is the parallel application 
of a quantum encoding algorithm [13] [5D] on fixed sized 
blocks, and a quantum computation on a fixed sized block 
can be efficiently simulated on a classical computer. 

Lemma 2. There exists a local encoding process £ , and 
decoding process, such that 
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for all times t and initial states \tp), i.e. simulation of 
the evolution of Hz is at least as hard as the simulation 
of the evolution due to Hi. 

Proof. In the proof of Lemma [T] we constructed an en- 
coding procedure £ on single spins that transformed the 
entire Hamiltonian Hi into a subspace of Hz- Thus, the 
action of H2 on any state confined within this subspace 
will be identical. One such example is the application of 
the same encoding procedure to any input state \ip). □ 

One might hope that the thermal state of H2 at inverse 
temperature (3 would accurately reproduce the properties 
of Hi at the same temperature, up to encoding and de- 
coding processes; we have designed H2 so that it can 
be written as a direct sum of Hi and other terms, all 
of which suffer energy penalties due to the added mag- 
netic field. This additional energy should correspond to 
a suppression of their importance in the thermal state, 
and hence not significantly affect the values of any ob- 
servables. Unfortunately, the vast number of these lev- 
els outweighs any potential saving, as we will now show. 
Consider an observable and its encoded counterpart 

h _ Tr(6 ie -^) - _ Tr(£(6i)e-^) 
[ 11 Tr(e-Pn) ' [ ' Tr(e-f> B ') ' 

We might expect some freedom in the way we define 
£{Oi). However, if we are to avoid obscuring the original 
observable, we must specify that £ {0\) = for any 
which is not a state encoded from the original space. 
Thus, Tr(£ (OOe - ^ 2 ) = T^Oie"^ 1 ), and hence 

(62) = Tr(e~^) 
(Oi) Tr(e-^)" 

This value is typically exponentially suppressed in N 
(this is readily proven for any J > max^ \ hi\), meaning 
that (O2) needs to be determined to exponential accuracy 
in order to rescale it to evaluate (Oi). Consider as an ex- 
ample a system with local magnetic field Hi = 
Each spin is encoded into 3 qubits to make it RI, so we 
write a local term of H 2 as (BZ) © ( Jl) where 1 is the 
6x6 identity matrix. 

(62) = ( 2cosh(E/3) \ N 
(Oi) V2co8h(B/3)+6e-0-V 

The denominator is clearly larger than the numerator 
due to the additional 6 levels, and the N th power effec- 
tively reduces the ratio to for large N and constant 
(3, J. While this does not constitute a general proof that 
observables cannot be determined, the fact that there 
must always be degeneracy present when translating Hi 
into H2 (since the quantum data cannot span different 
subsystems) implies that this is the case. 

Lemma 3. If Hi is constrained to a ID geometry and is 
TI, there also exists an instance H2 which is confined to 



a ID geometry, is TI and RI, and which recovers the pre- 
vious equivalences of ground state energy and dynamics. 

In this mapping, we will transform the period 1 TI of 
the d dimensional spins to a period 1 TI of qubits (i.e. the 
Hamiltonian acting on a set of qubits is the same for all 
translations along the chain) not a |~log 2 (eO] periodicity. 

Proof. As in Lemma [l] find d r > d and recode Hi, then 
apply the magnetic field. In the previous case, the Hamil- 
tonian terms ft-i...fc were transcribed to terms h' x rk . Due 
to the TI of Hi, all the terms h are the same, just acting 
on different blocks of qubits. However, now we will intro- 
duce a flag state \F), and a flag detection projector P F , 
both of F qubits. These RI flags will sit in front of each 
logical spin to herald the start of the logical spin. The 
state that the Hamiltonian acts on, and the projector 
P F , will be chosen such that P F has a non-zero outcome 
only if it aligns exactly with the qubits of the flag state. 
Introducing these states alters the Hamiltonian terms, 
such that the transcription takes the form 

^^■F+l...F+r,2F+r+1...2(F+r),...,(k~l)(F+r)+F+l...k(F+r)- 

Fig. [lja) depicts a suitable flag of the form 

pF _ p 2,0p2,lp2,0 p m,m/2 
"1...F ~ r l,3 r 2,4 r 5,6 r 7...6+rai 

for some even integer m, to be selected momentarily. It 
is also necessary to introduce a further magnetic field 

N(F+r)-l 

J' E T ^- P l...F®P r FU...F + r)T l 

i-0 

where J' = 2fc'maxj|/^| so that states that are not 
of the correct form are penalised. The Hamiltonian is 
k' = k(F + r)-body acting on (F + r)N qubits. With 
periodic boundary conditions, the eigenstates of H2 are 
(3(m+ 1)) N (F + r) times more degenerate than Hi, cor- 
responding to the different possible alignments of the 
ground state, and the different choices of state detected 
by P 2 ' 1 , of which there are 3 per flag state, and the 
rn+1 different choices for the p m >™/ 2 term. It remains to 
prove that C = T~ n Pf F T n , a term which will appear in 
the Hamiltonian, has no support on AB, the flag states 
from two neighboring encoded spins in the (ground) state 
of the system. These flag states constitute any of the 
rank(P f ) orthogonal states that satisfy P F \F) — \F). 
By taking A = P[ F and B = T- {F+r ^P F „ F T F+r , we 
require ABC = for n = 1 . . . r + F — 1. Consider, first, 
the case of n = 1. Here the term in C overlaps with 

P2 4 in A, and thus ABC = P24P2 I = 0- Similarly, the 
P5 ' 6 term in C hits p m >™ 1 / 2 m a for n — 2 . . . m, and the 
P2 '4 also hits over the range n — 6 . . . m + 3, resulting 
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FIG. 1: The flag projector P . There is a duality between 
the projectors in the Hamiltonian and the ground state of 
the Hamiltonian - one selects any of the states that the flag 
projector projects onto, (a) A general construction that works 
for any values of r and j for sufficiently large m. (b) Special 
case construction to give F ~ 0(1) for j m r/2. 

in overlap, Eqn. pi). There is a symmetry between the 
overlaps of C and A at offset n and those of B and C at 
offset r + F — n. By imposing that m = max(r — 1,5), 
these offsets cover the entire range. 

Smaller values of m are often possible depending on 
the Afj subsystem that we encode into. For example, 
if r/2 + j + 1 < m < r, then for offsets n — m to r, 
the p m > m / 2 in C entirely overlaps the encoded spin of 
A, which is in the j subspace, and the projection is 0, 
Eqn. (JT|). When taken in conjunction with the symmet- 
ric version on B and C, then all offsets n give projection 
provided r/2 > j + 6, and we select m = r/2 + j + 1. In 
the limit of large r, with j — \/r/2, the most efficient en- 
coding into a decoherence-free subsystem, one finds that 
the cost of a TR1 Hamiltonian scheme compared to a Tl 
scheme is an increase in the range of the Hamiltonian 
terms of approximately 50%. For the smallest possible 
values of r = 3,4, P F can be altered, replacing P^P^'l 
with P13P2 4 and setting m — r+1, which yields a minor 



improvement. For large values of j (of order r/2), one 
can find superior choices for the flag state, of constant 
size. However, there is little interest in elucidating this 
construction, and is relegated to the example in Fig.Qb). 

ft is immediate from the preceding arguments that 
there exists an encoding procedure £ to satisfy the TI 
variants of Lemmas [T] and [2] □ 

Upon application of this result to the Hamiltonian of 
|21j . which required TI nearest-neighbor interactions of 
10-level systems to implement universal quantum com- 
putation through dynamics, the result is a 38-body TRI 
local ID Hamiltonian (r = 7,j — \,m — 6). This is 
compared to the 106-body interactions required in [§]• 

Conclusions: Imposing RI on a Hamiltonian does not 
affect the difficulty of finding the ground state energy - 
for any Hamiltonian, one can always write down a RI 
version that accurately reproduces all salient properties, 
including, in particular, TI. Hence it is not expected that 
TRI Hamiltonians which are local on a ID lattice can 
be efficiently solved [11]. The same can be said for the 
iV-representability problem as applied to a fc'-qubit den- 
sity matrix p [TTJ [22] , which should be compatible with 
a TRI pure state \tp) of N qubits. The results of [23] 
can also be applied in order to discuss the universality of 
Hamiltonians for Adiabatic Quantum Computation. Un- 
fortunately, none of the results apply to the calculation of 
expectation values of thermal states at any finite temper- 
ature. It remains an interesting open question whether 
alternative formulations can prove the persistence of the 
complexity of calculation for finite temperatures. 

Acknowledgments: N. Schuch is thanked for a critical 
reading of the manuscript . This work is supported by the 
DFG Cluster of Excellence Munich-Centre for Advanced 
Photonics (MAP) and Clare College, Cambridge. 



[1] B. A. Chase and A. J. Landahl, |arXiv:0802.1207| (2008). 

[2] D. Aharonov et al. Proceedings of the 48th IEEE Sym- 
posium on Foundations of Computer Science, FOCS 07 
(IEEE, Providence, RI, 2007), p. 373. 

[3] A. Y. Kitaev, A. H. Shen, and M. N. Vyalyi, Classical and 
quantum computation, Graduate studies in mathematics 
(AMS, Providence, Rhodes Island, 2002). 

[4] R. Oliveira and B. M. Terhal, Quant. Inf. Comp. 8, 0900 
(2008). 

[5] J. Kempe, A. Kitaev, and O. Regev, SIAM Journal of 

Computing 35, 1070 (2006). 
[6] J. D. Biamonte and P. J. Love, Phys. Rev. A 78, 012352 

(2008). 

[7] M. B. Hastings, Phys. Rev. B 73, 085115 (2006). 

[8] K. G. H. Vollbrecht and J. I. Cirac, Phys. Rev. Lett. 100, 

010501 (2008). 
[9] A. Kay, Phys. Rev. A 78 012346 (2008). 
[10] A. Kay, Phys. Rev. A 76, 030307(R) (2007). 
[11] D. Gottesman and S. Irani, arXiv: 0905.2419 (2009). 
[12] P. Zanardi, and M. Rasetti, Phys. Rev. Lett. 79, 3306 
(1997). 



[13] P. Zanardi, Phys. Rev. Lett. 87, 077901 (2001). 

[14] E. Knill, R. Laflamme, and L. Viola, Phys. Rev. Lett. 
84, 2525 (2000). 

[15] S. D. Bartlett, T. Rudolph, and R. W. Spekkens, Phys. 
Rev. A 70, 032307 (2004). 

[16] S. D. Bartlett, T. Rudolph and R. W. Spekkens, Rev. 
Mod. Phys 79, 555 (2007). 

[17] M. S. Byrd, Phys. Rev. A 73, 032330 (2006). 

[18] S. J. van Enk, Phys. Rev. A 73, 042306 (2006). 

[19] D. Bacon, I. L. Chuang, and A. W. Harrow, Phys. Rev. 
Lett. 97, 170502 (2006). 

[20] D. Bacon, I. L. Chuang, and A. W. Harrow, Proceed- 
ings of the eighteenth annual ACM-SIAM symposium on 
Discrete algorithms (Society for Industrial and Applied 
Mathematics, Philadelphia, 2007), p. 1235. 

[21] D. Nagaj and P. Wocjan, Phys. Rev. A 78, 032311 (2008). 

[22] Y. K. Liu, M. Christandl, and F. Verstraete, Phys. Rev. 
Lett. 98, 110503 (2007). 

[23] D. Aharonov et al. SIAM Journal of Computing, 37, 166 
(2007). 



